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ABSTRACT 

We investigate properties of the interstellar medium (ISM) in galaxies hosting long-duration 7-ray 
bursts (GRBs) from an analysis of atomic species (Mg", Fe") and excited fine-structure levels of ions 
(e.g. Si+). Our analysis is guided primarily by echelle observations of GRB 050730 and GRB 051111. 
These sightlines exhibit fine-structure transitions of O", Si+, and Fe+ gas that have not yet been detected 
in intervening quasar absorption line systems. Our results indicate that the gas with large Mg I equivalent 
width (e.g. GRB 051111) must occur at distances > 50 pc from GRB afterglows to avoid photoionization. 
We examine the mechanisms for fine-structure excitation and find two processes can contribute: (1) 
indirect UV pumping by the GRB afterglow provided a far-UV intensity in excess of 10^ times the 
Galactic radiation field; and (2) coUisional excitation in gas with electron density > lO^cm"^. The 
observed abundances of excited ions are well explained by UV pumping with the gas at r a few hundred 
pc from the afterglow for GRB 051111 and r < 100 pc for GRB 050730, without invoking extreme gas 
density and temperature in the ISM. We show that UV pumping alone provides a simple explanation for 
all reported detections of excited ions in GRB afterglow spectra. The presence of strong fine-structure 
transitions therefore may offer little constraint for the gas density or temperature. We discuss additional 
implications of UV pumping including its impact on chemical abundance measurements, new prospects 
for observing line-strength variability, and future prospects for studying the gas density and temperature. 
Finally, we list a series of criteria that can distinguish between the mechanisms of UV pumping and 
coUisional excitation. 



Subject headings: gamma-ray bursts 
1. INTRODUCTION 

Optical afterglows of long-duration 7-ray bursts (GRBs) 
are almost certainly signposts of extreme star-forming re- 
gions at high redshift, because these bursts originate in 
the catastrophic death of massive stars (e.g. Woosley 1993; 
Paczyhski 1998; Bloom et al. 2002; Stanek et al. 2003). 
A detailed study of metal absorption features identified in 
the circumburst environment can, in principle, yield strong 
constraints on the progenitor (e.g. Perna & Loeb 1998; 
Ramirez-Ruiz et al. 2005; van Marie, Langer & Garcia- 
Segura 2005) and offer important insights for understand- 
ing mass-loss and chemical feedback during the final evo- 
lution stages of massive stars (Wijers 2001; Vink 2005). 
The benefit of GRBs, which presumably are embedded in 
environments of Wolf-Rayet stars, is that their optical af- 
terglows provide extremely bright lighthouses at moderate 
to high redshift, allowing redshifted UV lines to be (albeit 
transiently) observable with high-resolution spectroscopy. 

Past studies based on low-resolution afterglow spectra 
have been limited to a few diagnostic measurements such 
as the H I column density, the gas metallicity, and the 
dust-to-gas ratio (e.g. Savaglio, Fall & Fiore 2003; Vreeswijk 
et al. 2004). In contrast, high-resolution spectroscopy of 
GRB afterglows has uncovered detailed kinematic signa- 
tures, population ratios of excited ions, and chemical com- 
positions of the interstellar medium (ISM) and, perhaps, 
the circumstellar medium (CSM) of the progenitor star of 
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the GRB (Vreeswijk et al. 2004; Chen et al. 2005). These 
quantities, in principle, permit direct comparisons with 
simulations describing the evolution history of GRB pro- 
genitors (e.g. van Marie, Langer & Garcia-Segura 2005). 

An emerging feature of GRB progenitor environments is 
the presence of strong fine-structure transitions from ex- 
cited states of C+, Si+, 0°, and Fe+ (Chen et al. 2005). 
Together with their resonance transitions, these fine-structure 
lines can in principle provide robust constraints on the 
temperature and density of the g well as the ambient 
radiation field, independent of the ionization fraction, rela- 
tive abundances, or metallicity of the gas (Bahcall & Wolf 
1968). In particular, detections of Fe II fine structure tran- 
sitions have only been reported in rare places such as broad 
absorption-line (BAL) quasars (Hall et al. 2002), r] Cari- 
nae (Gull et al. 2005; Nielsen, Gull & Vieira Kober 2005) 
and the circumstellar disk of (3 Pictoris (Lagrange-Henri, 
Vidal-Madjar & Ferlet 1988). The presence of strong Fe II 
fine-structure lines therefore suggests extreme gas density 
and temperature in the GRB progenitor environment. In 
the case of circumburst medium, however, the greatly en- 
hanced radiation field due to the GRB afterglow may also 
contribute significantly to the formation of these excited 
ions. 

In this paper we present a comprehensive study of the 
physical properties of the ISM in GRB host galaxies, us- 
ing absorption lines identified in the echelle spectra of two 
GRB afterglows, GRB 051111 and GRB 050730 as case 
studies. We first consider the enhanced radiation field in 
the circumburst environment due to the progenitor star 
as well as the optical afterglow, and constrain the dis- 
tance of the observed neutral gas cloud based on the pres- 
ence/absence of various atomic lines such as Mg I and Fe I. 
The dominant species of these elements is in the singly ion- 
ized state, because their first ionization potential is less 
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Fig. 1. — Velocity profiles for transitions arising from gas in the host galaxy of GRB 051111. The vertical dashed line at d = corresponds 
to z = 1.5495 and the dotted lines indicate blends with coincident transitions. 
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than 1 Ryd. The presence of strong atomic Unas there- 
fore places a lower limit to the distance of the cloud from 
the afterglow. Next, we investigate the excitation mecha- 
nism of the absorbing gas, again taking into account the 
presence of an intensified radiation field known from after- 
glow light-curve observations. We compare the predicted 
population ratios with observations of different ions under 
photon pumping and coUisional excitation scenarios. Al- 
though the analysis is based on echelle data obtained for 
two GRB hosts, we note that the results are applicable to 
the majority of long-duration GRB's discovered to date. 
Finally, we discuss additional framework for studying the 
ISM gas of future events. 

The paper is organized as follows. The observations and 
analysis of echelle spectra for GRB 051111 and GRB 050730 
are briefly summarized in § 2. In § 3 we describe the 
radiation field produced by GRB afterglows. In § 4 we 
analyze the non-dominant atomic transitions observed in 
optical afterglow spectra, while § 5 investigates the mecha- 
nisms for populating fine-structure states in the dominant 
ions. In § 6 we consider tests to distinguish between UV 
pumping and coUisional excitation as the principal mech- 
anism for population of fine-structure levels. Finally, § 7 
presents a discussion of the results for the specific cases of 
GRB 051111 and GRB 050730 and, also, generic results 
for observations of long-duration GRBs. Throughout the 
paper we assume a cosmology with il\ = 0.7, flm = 0.3 
and Ho = 70kms"iMpc"^ 
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Fig. 2. — Velocity profiles for Fe II transitions arising from gas in 
the host galaxy of GRB 050730. The vertical dashed line at = 
corresponds to z = 3.96855 and the dotted lines indicate blends with 
coincident transitions. 



2. OBSERVATIONS AND COLUMN DENSITIES 

Echelle spectroscopic observations of GRB 051111 were 
obtained using the High Resolution Echelle Spectrometer 
(HIRES; Vogt et al. 1994) on the Keck I telescope. The 
observations consisted of three exposures of 1800 s du- 
ration each and were initiated by the Keck Observatory 
staff at UT 07:02:57.92 on 11 November 2005. The spec- 
tra^ covered a wavelength range from A = 4165 - 8720A 
with a spectral FWHM resolution of « 5 km s~^. A sum- 
mary of the data reduction is presented in a separate paper 
(Prochaska et al. 2006, in preparation). The final coad- 
ded spectrum was normalized to a unit continuum and 
has SNR « 15 per 1.3 km s""'^ pixel. Figure 1 presents the 
subset of transitions related to the GRB for the analysis 
performed in this paper. We (Hill et al. 2005) reported 
the discovery of the redshift z = 1.549 for GRB 051111 
and preliminary analysis of the results Prochaska (2005). 
Penprase et al. (2005) present an independent analysis of 
the same data. 

Echelle spectroscopic observations of GRB 050730 were 
obtained using the MIKE echelle spectrograph (Bernstein 
et al. 2003) on the Magellan Clay telescope. The data 
cover observed wavelengths A = 3400-9000A at a FWHM ? 
11km s^^ dispersion with a typical SNR=12 per resolu- 
tion element. A detailed description of the observation and 
data reduction is presented in Chen et al. (2005). In ad- 
dition to the absorption features presented in Chen et al., 
here we also present identifications of Fe II transitions from 
excited fine structure states, 3c?^4s ^-D ,7=5/2, 3/2, 1/2 (Fig- 
ure 2). Column density measurements of the two GRB 
sightlines are presented in Tables 1 and 2. The atomic 
data are taken from Morton (1991); Bergeson & Lawler 
(1993); Bergeson, MuUman & Lawler (1994, 1996); Tripp, 
Lu & Savage (1996); Raassen & Uylings (1998); Morton 
(2003). 

Column densities of the unblended transitions were mea- 
sured using the apparent optical depth method (AODM; 
Savage & Sembach 1991; Jenkins 1996) over a velocity in- 
terval ±40 km s^i centered at z = 3.96855 for GRB 050730 
and at 2 = 1.54948 for GRB 051111. In Tables 1 & 2, we 
list the column densities of individual transitions, the Icr 
statistical uncertainty, and the weighted mean. In addi- 
tion, we report conservative lower limits for saturated tran- 
sitions and 3f7 upper limits for non-detections. No system- 
atic trend between column density and oscillator strength 
is identified in the measurements, indicating that the ab- 
sorption lines are well resolved and that the column density 
measurements are robust. We note, however, that uncer- 
tainties in the continuum level are large for features with 
rest-frame absorption equivalent width < 10 mA. In spite 
of a formal statistical error < 0.01 dex in several cases, we 
adopt a minimum la error of 5% in subsequent analyses. 

3. RADIATION FIELD PRODUCED BY THE GRB 
AFTERGLOW 

The intimate connection between long-duration GRBs 
and the end stage of the lifetime of massive stars implies 
that the gaseous clouds in the vicinity of the progenitor 
stars have been blasted with intense ultraviolet radiation 
fields, especially in the early afterglow phase. The radi- 
ation field near the GRB, significantly enhanced with re- 

*The data are publicly available at http;//www. graasp.org 
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Table 1 

IONIC COLUMN DENSITIES FOR GRB 051111 



Ion 




cm~^ 


A 

(A) 


/ 


logN 


log ^adopt 


Mgl 


















0.00 


1747.794 


0.009080 


14.68 ± 0.04 


14.68 ± 0.04 






0.00 


1827.935 


0.023900 


> 14.63 








0.00 


2026.477 


0.112000 


> 14.31 




Aim 


















0.00 


1862.790 


0.268000 


13.51 ±0.02 


13.51 ±0.02 


Si II 
















1/2 


0.00 


1808.013 


0.002186 


> 16.07 


> 16.07 




1/2 


0.00 


2335.123 


0.000004 


< 17.27 






3/2 


287.24 


1816.928 


0.001660 


15.00 ±0.06 


15.00 ±0.06 


Fel 
















0.00 


2484.021 


0.557000 


< 11.75 


< 11.75 


Fell 
















9/2 


0.00 


2249.877 


0.001821 


15.20 ±0.02 


15.24 ±0.01 




9/2 


0.00 


2260.780 


0.002440 


15.28 ± 0.02 






7/2 


384.79 


2333.516 


0.069170 


13.97 ± 0.01 


13.98 ±0.01 




7/2 


384.79 


2365.552 


0.049500 


13.99 ± 0.01 






7/2 


384.79 


2383.788 


0.005175 


14.05 ± 0.09 






7/2 


384.79 


2389.358 


0.082500 


14.01 ±0.02 






7/2 


384.79 


2626.451 


0.044100 


< 14.00 






5/2 


667.68 


2328.111 


0.035550 


13.75 ±0.02 


13.76 ± 0.01 




5/2 


667.68 


2607.866 


0.118000 


13.76 ± 0.02 






5/2 


667.68 


2618.399 


0.050500 


13.76 ± 0.02 






3/2 


862.61 


2338.725 


0.092500 


13.54 ±0.02 


13.55 ±0.01 




3/2 


862.61 


2359.828 


0.057270 


13.59 ±0.02 






3/2 


862.61 


2621.191 


0.003920 


< 13.86 






1/2 


977.05 


2345.001 


0.157500 


13.17 ±0.03 


13.25 ± 0.01 




1/2 


977.05 


2414.045 


0.175500 


13.25 ±0.02 






1/2 


977.05 


2622.452 


0.056000 


13.28 ± 0.04 






1/2 


977.05 


2629.078 


0.173000 


13.28 ±0.02 




Zn II 
















0.00 


2026.136 


0.489000 


> 13.60 


> 13.71 






0.00 


2062.664 


0.256000 


> 13.71 





"J value for ions with excited states. Ei^^ is the energy above the ground 

state. 

'"Although the statistical uncertainty is as low as 0.01 dex in several cases, 
we adopt a minimum uncertainty of 5% in the analysis. 



spect to the ambient radiation intensity of the host galaxy, 
modifies the ionization and excitation states of the gas. It 
is therefore critical to consider first the radiation field pro- 
duced by the GRB afterglow before then analyzing the 
absorption-line profiles. 

Wc adopt the afterglow of GRB 051111 as an example. 
Butler et al. (2006) have analyzed the light curve data 
of this event and characterize the flux of the afterglow at 
Earth for tohs > 50 s as 

Ff« = 1.4 X 10-25 (^-^\ "7 ''ergcm-^s-iRz 

(1) 

where Aobs and iobs are the observed wavelength and time 
at Earth. The fading of the afterglow is best-fit with 
at = 0.87. The spectral index was estimated to be /3 = 0.6 
from a fit to the VRIJHK bands, i.e. < 6cV at the rest 
frame of the GRB host. If we assume that this power-law 
extends to higher energies, then Equation 1 describes the 
unattenuated flux to energies of several Ryd as a function 
of time. At some level this radiation field is attenuated by 
(1) dust and gas in the host galaxy: Lyman limit absorp- 
tion at hf > IRyd, (2) dust extinction at all wavelengths, 
and (3) line opacity by heavy elements and hydrogen Ly- 



man series (concentrated at energies hu > 6eV). All of 
these sources of opacity arc expected to contribute and we 
will discuss them when relevant. 

We consider three energy intervals important to UV/optical 
absorption line spectroscopy: (1) a lA interval centered at 
hiy = 0.048 cV (385cm-i);' (2) 8eV < hiy < 13.6 eV; and 
(3) IRyd < hiy <2Ryd. The first is set by the energy 
necessary to populate the first excited state J = 7/2 of 
Fe"*" from its ground state J = 9/2 via a magnetic dipole 
transition. The lA interval is chosen to roughly coincide 
Vth the width of the line-profile which is set by the gas 
kinematics along the GRB sightline. The second inter- 
val is a measure of the far-UV radiation field. To provide 
comparison with the Milky Way, we adopt Habing's con- 
stant Go = 1.6 X 10""^ ergcm-2 s"^, which roughly cor- 
responds to the Galactic far-UV intensity (Habing 1968). 
More recent estimates place the Galactic radiation field 
at « 1.7Go (Gondhalckar, PhiUips & Wilson 1980). The 
third energy interval characterizes the radiation field re- 
sponsible for ionizing hydrogen atoms. 

In Figure 3 we plot the number of photons per cm^ emit- 
ted by the afterglow of GRB 051111 from t = tohs/il + z) = 
19s to 1506s versus the radial distance r from the GRB in 
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Table 2 

IONIC COLUMN DENSITIES FOR GRB 050730 



Ion 


J" 




A 

(A) 


/ 


logN 


log Nadopt 


CI 
















1/2 


0.00 


1277.245 


0.096650 


< 13.37 


< 13.37 




1/2 


0.00 


1656.928 


0.140500 


< 13.47 




OI 














1/2 


0.00 


1302.168 


0.048870 


> 15.27 


> 15.27 




1/2 


0.00 


1355.598 


0.000001 


< 18.09 






2/2 


158.26 


1304.858 


0.048770 


> 14.94 


> 14.94 




4/2 


226.98 


1306.029 


0.048730 


> 14.61 


> 14.61 


Si II 














1/2 


0.00 


1260.422 


1.007000 


> 14.08 


> 14.69 




1/2 


0.00 


1304.370 


0.094000 


> 14.69 






1/2 


0.00 


1526.707 


0.127000 


> 14.55 






1/2 


0.00 


1808.013 


0.002186 


< 15.59 






3/2 


287.24 


1264.738 


0.903400 


> 14.03 


< 16.03 




3/2 


287.24 


1265.002 


0.100400 


> 14.56 






3/2 


287.24 


1309.276 


0.146800 


> 14.28 






3/2 


287.24 


1817.451 


0.000129 


< 16.03 




Fell 














9/2 


0.00 


1608.451 


0.058000 


> 14.72 


14.98 ±0.12 




9/2 


0.00 


1611.200 


0.001360 


14.98 ±0.12 






7/2 


384.79 


1618.468 


0.021400 


14.26 ±0.06 


14.29 ±0.04 




7/2 


384.79 


1621.686 


0.038100 


14.34 ±0.07 






5/2 


667.68 


1629.160 


0.036700 


14.12 ±0.07 


14.12 ±0.07 




3/2 


862.61 


1634.350 


0.020500 


13.95 ±0.11 


13.92 ±0.06 




3/2 


862.61 


1636.331 


0.040500 


13.91 ±0.07 






1/2 


977.05 


1639.401 


0.057900 


13.65 ±0.07 


13.65 ±0.07 



"^J value for ions with excited states. -E;om is the energy above the ground 
state. 

'Although the statistical uncertainty is as low as 0.01 dex in several cases, 
we adopt a minimum uncertainty of 5% in the analysis. 



the rest-frame of the GRB host. For the E — 385 cm"-'^ 
photons, we assume a lA interval. Similar to previous 
work (e.g. Perna & Loeb 1998; Draine & Hao 2002), we find 
that in the absence of intervening gas and dust the GRB 
afterglow would initiate a burning front that would ionize 
an H I column of several 10^^ cm~^ atoms at r = 1 pc. 
As impressive, one notes that the far-UV intensity is over 
10^ times larger than the ambient Milky Way field even 
at r = 1 kpc. Despite its transient nature, the afterglow 
emits more UV photons in a few minutes than the entire 
galaxy in one year (assuming a star formation rate SFR 
= 20 Mq yr~i = 1.6 x lO^^ergs"^ Hz"!; e.g. Kennicutt 
1998). 

The estimated strength of the radiation field due to this 
afterglow has three important implications, which are ap- 
plicable to many other GRBs. First, if neutral gas were 
present at r < 1 pc of the progenitor prior to the burst 
(which is unlikely assuming significant UV luminosity from 
the GRB progenitor), the implied mean hydrogen density 
< Uh > would have to exceed lO^cm"'^ for any neutral 
gas to survive the burst. At r > 10 pc, however, significant 
absorption is still expected from neutral gas if < nn > is 
greater than 10 cm"'^. Second, neutral atoms with ioniza- 
tion potential less than 1 Ryd (present within the spec- 
trum of GRB 051111) are very unlikely to survive the first 
few hundred seconds (§ 4). Third, the intensity of the 
radiation field at even 1 kpc is large enough that photon 
pumping of fine-structure states cannot be ignored. In the 
following sections, we will consider photon pumping and 




1 10 100 1000 

Distance from GRB (pc) 



Fig. 3. — Integrated, unattenuated flux of photons emitted by the 
afterglow of GRB 051111 from t = 19.6 to 1506s with energies in 
three intervals: (1) Hydrogen ionizing radiation (IRyd < hu <2Ryd; 
red curve); (2) far-UV radiation (8eV < hu <13.6eV; black curve); 
(3) photons which would excite Fe"*" from the J = 9/2 to J = 7/2 
level (a lA interval centered at hu = 0.048eV; blue curve). The 
radial dependence is simple dimming. The right axis presents the 
time-averaged intensity of the far-UV radiation field Ggrb relative 
to Habing's constant Go. 
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discuss constraints on the distance to the progenitor stars, 
the gas density and temperature of the observed gas. 

4. GASEOUS PROPERTIES AS CONSTRAINED BY ATOMIC 

SPECIES 

By definition, neutral hydrogen gas clouds are opaque 
to radiation with hiy > 1 Ryd and transparent to lower 
energy photons (aside from opacity by the Lyman series). 
Because the first ionization potential of most elements is 
less than 1 Ryd, the dominant species of these elements 
is singly ionized gas, e.g. Fe+, Si+, and Mg+. However, 
the atomic phase is detected in trace quantities in various 
gaseous phases of the ISM of nearby galaxies (e.g. Frisch 
et al. 1990; Welty, Hobbs & Morton 2003). One of the 
most common examples is Mg*^. Surveys of Mg I lines 
in the Milky Way have revealed this atom in both cold 
and warm phases of the neutral ISM (e.g. Frisch, York & 
Fowler 1987; Frisch et al. 1990). In addition, atomic Mg 
is also often detected in quasar absorption line systems 
exhibiting strong Mg II absorption (e.g. Stcidcl & Sargent 
1992; Churchill et al. 2000). This includes damped Lya 
systems at all redshifts (Rao & Turnshek 2000; Prochaska 
& Wolfe 2002). The frequent presence of Mg° is driven by 
a relatively low photoionization cross-section and modest 
recombination coefficients. At low temperature radiative 
recombination dominates while at T > 5000K dielectronic 
recombination becomes significant, allowing Mg I absorp- 
tion to be observed from even partially ionized regions 
(York & Kinahan 1979). 

In this section, wc examine the presence/ absence of Mg" 
(IP=7.64eV) and Fe° (IP=7.90eV) along the GRB 051111 
sightline. Although the results are directed to the gas in 
this GRB host galaxy, neutral species have been reported 
in other GRB afterglow spectra (e.g. Barth et al. 2003; 
Foley et al. 2005) and our conclusions are qualitatively ap- 
plicable to those cases. One can also generalize the results 
to additional atomic species, e.g., Si'^ and C^. 

Column density mcasiircments for the atomic species 
shown in Figure 1 are presented in Table 1. We detect 
several transitions of Mg I but no Fe I transition. The 
Mg° column density for the host galaxy of GRB 051111 
may be the largest ever recorded outside the Milky Way. 
Based on the kinematics of the Mg I profile and its large 
column density, we argue that this atomic gas is co-spatial 
with the majority of low-ion absorption observed along the 
sightline. As such, we will treat the Mg°, Mg+, Fc", and 
Fe"*" gas to arise from a single cloud within the ISM of 
the galaxy. While Mg° atoms are also likely to reside in 
the halos of galaxies (e.g. Charlton & Churchill 1998), the 
expected column density is small and should show signifi- 
cantly diflferent kinematics. 

4.1. Distance of Neutral Gas to the GRB 

Even if physical conditions in the ISM of the GRB host 
galaxy support a significant Mg" column density, the en- 
hanced far-UV radiation from the GRB afterglow can burn 
away the gas even at large distance from the burst. Ob- 
servations of Mg°, therefore, can place a lower limit to the 
distance of the neutral gas to the burst. 

As an estimate to this lower limit, we calculate the dis- 
tance rmin at which 99.99% of the Mg" atoms would be 
ionized by the GRB afterglow. At this level of photoion- 
ization, one would not expect to detect even the strong 



Mg I 2852 transition. For an optically thin slab, rmin cor- 
responds to the distance at which the gas experiences an 
integrated photon surface density = ln(10^)/crpft. For 
a total number of ionizing photons 0, the distance is 



100 X 



1.1 X 10^3 



pc 



(2) 



We adopt the functional form given by Verner et al. (1996) 
for the photoionization cross-section of Mg" , (T^^ (E). This 
function peaks at the ionization potential 
a^/{E = 7.64eV) = 1.7 x lO'^^ cm^ and declines sharply 

with increasing energy. We assume a value ct^^ = 10~^^ cm^ 

(corresponding to £" w 8.5eV) as a conservative estimate 
of the cross-section. For GRB 051111, wc calculate a time- 
averaged luminosity density at hv = 8.5cV, 
< >j= 6.4 x 10^^ ergs~^ cm~2 Hz~i over t = 20 - 
1000 s, leading to a mean total photon number = 2 x 10^*^ 

at 8.5 eV. Substituting cr^^ for aph in Equation (2) yields 

rmin = 46 pc for the neutral gas observed in the GRB 051111^. 
The presence of strong Mg I absorption transitions at 
t = 1500s after the burst therefore suggests that the neu- 
tral gas is located at least 50 pc away from the afterglow. 



Table 3 

Constraints on Circumburst Distances of Observed 

Neutral Gas 



GRB 


2 


a 


13 


Rcf 


log Li 


Mgl 


excite 












(cgs) 


(pc) 


(pc) 


010222 


1.477 


0.80 


0.89 


1 


31.39 


40 


190 


020813 


1.254 


0.85 


0.92 


2 


31.09 


30 


140 


021004 


2.328 


1.05 


1.05 


3 


32.21 


140 


620 


030323 


3.372 


1.56 


0.89 


4 


32.85 


540 


2330 


030329 


0.169 


1.10 


1.00 


5 


31.38 


60 


250 


050408 


1.236 


0.79 


1.30 


6 


29.93 


10 


40 


050730 


3.969 


0.30 


1.80 


7 


32.16 


70 


340 


050820 


2.615 


0.95 


1.00 


8 


31.97 


100 


430 


051111 


1.549 


0.87 


0.60 


9 


31.32 


40 


180 


060206 


4.048 


1.01 


0.51 


10 


32.41 


170 


730 



"Specific luminosity of the GRB afterglow at t = 1000s for i/ = 8eV//i. 

''Radius from the GRB afterglow where 99.99% of Mg" atoms 
would be ionized. 

"^Approximate radius where a column of 10^* cm~^ Si+ ions would 
experience one excitation per minute per ion. 

Note. — We caution that the values listed in this Table should 
be considered rough estimates of the light curves. 
The flux of the GRB observed at Earth is given by 
F„ = (L„/47r(i2/[l + 2])(t/1000)«(i//1.7 x IQi^Hz)''. 
We have not corrected for the fluxes for dust extinction 
in the host galaxy. 

References. — 1; Mirabal et al. (2002); 2; Barth et al. (2003); 
3: Mirabal et al. (2003); 4: Vreeswijk et al. (2004); 
5: Lipkin et al. (2004); 6: Foley et al. (2005); 
7: Chen et al. (2005); 8: Cenko & Fox (2005); 
9: Butler et al. (2006); 10: Perley et al. (2006) 



A stronger constraint can be placed on the minimum 
distance from the GRB by evaluating time variations in the 

^Note that one can ignore recombination provided that the den- 
sity is not extraordinary, i.e. Tie < lO^cm"^. 
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Mg° column density (Mirabal et al. 2002). In this case, the 
result is independent of the total Mg° column density prior 
to the GRB. For GRB 051111, we set an upper limit to the 
decline in A^(Mg*') over the central 40km s~^ interval to 
be less than 25% during a Aiobs = 1800 s exposure. The 
number of photons emitted over the time interval tobs = 
(3840 + 900) s to tobs = (3840 + 2700) s is = 2.3 x lO". 
Therefore, we find a lower limit to the separation based on 
the absence of line-strength variability to be 



Mg" 
'ph 



47r \n{Ni/Nf) 



> 



I pc 



(3) 



for Nf > 0.75Ni. This is in good agreement with our 
estimate for r^i^ from above. Note that if we were able 
to place an upper limit to the variability of 5% then r^^^ 
would double. 

We present estimates of Tmin for several other GRB as- 
suming Equation 2, together with the estimated parame- 
terization of the light curves shown in Table 3. The results 
for GRB 020813 and GRB 050408 are notable because 
those sightlines also exhibit significant Mg° gas. Our anal- 
ysis demonstrates that the conditions for GRB 051111 are 
not atypical. A generic result of our analysis is that Mg" 
absorption identified in GRB spectra must correspond to 
gas at distances > 100 pc from the afterglow. 



The corresponding lower limit of rie x (G/Gq) "'^ as a func- 
tion of T is shown in the solid, blue curve in Figure 4. 



0.100 



o 
o 



0.010 



0.001 




10000 



Fig. 4. — Constraints on the ratio of electron density ris to the 
ambient far-UV intensity (G/Gq) as a function of the gas tempera- 
ture. The black curve set an upper limit based on the non-detection 
of atomic Fe. The blue curve set a lower limit based on the detection 
of atomic Mg and a conservative upper limit to the column density 
of Mg+. 



4.2. Electron Density 

In the previous section, we have argued that the very 
presence of Mg° gas or lack of time variability in its column 
density indicate the gas is located at large distance from 
the GRB afterglow. In this case, one can examine the gas 
under the assumptions of steady-state balance and pho- 
toionization equilibrium. Considering together the pres- 
ence (or absence) of atomic species like Mg° and Fe° in 
a neutral hydrogen gas leads to constraints on the ratio 
of the electron density to the intensity of ambient far-UV 
radiation. For the ionization balance of Mg" and Mg"*", we 
express photoionization equilibrium in terms of the column 
density ratios (Appendix A), 

^ jV(MgO) r(Mg") 

where r(Mg°) = Const, x aph{Mg°) x G/Gq represents 

the photoionization rate of Mg° and a^^^ (T) is the total 
recombination coefficient dependent upon the gas temper- 
ature T. Equation (4) demonstrates that observations of 
the relative ionization fraction for a single heavy element 
constrain n^ x (G/Gq)~^ at a given temperature. 

For GRB 051111, the Mg II doublet is highly saturated 
in the echelle data and only a lower limit can be placed. 
We can estimate, however, a conservative upper limit to 
A''(Mg"'") from the observed column density of Zn+. First, 
we assume relative solar abundances for Zn and Mg. Next, 
we increment the observed log 7V(Zn+) by 1 dex to account 
for both line-saturation and a nucleosynthetic enhance- 
ment relative to solar. This treatment presumes the Mg*' 
gas is co-spatial with the Mg+ ions inferred from Zn+, 
which is supported by the kinematic characteristics of the 
various line profiles (Figure 1). Finally, we derive a conser- 
vative lower limit to the ratio 7V(Mg°)/iV(Mg+) > lO"^-*. 



We repeat this analysis for the observed ratio Fe°/Fe"'". 
At T < 8000K, radiative recombination dominates and we 
have 

_ A^(Fe°) r(Fe") 

7V(Fe+)aF<=+(T) ' ^ ' 

The absence of Fe*' places an upper limit on the column 
density ratio iV(Fe*')/iV(Fe'*') < lO^'^, which in turn leads 
to an upper limit of rie x (G/Gq)^^ as shown in the solid, 
black curve in Figure 4. 

Figure 4 allows us to roughly estimate the electron den- 
sity of the gas for a known ambient radiation field G. If 
the GRB progenitor was located within a star forming 
region, then the ambient far-UV radiation field could be 
large. Observations of known GRB host galaxies at z > 1 
show a mean star formation rate integrated throughout 
the host ISM of SFR- WM^yr'^ (Christenscn, Hjorth 
& Gorosabel 2004; Le Floc'h et al. 2006). Even if all 
the star formation were concentrated in a single cloud, the 
mean intensity at Ikpc is G/Gq ~ 100. At T < lOOOOK, 
we therefore estimate = 1 — 10cm~^. Assuming an 
ionization fraction of a; = nH+/nH = 10"^, we estimate 
the hydrogen volume density nn — 100 — 1000 cm~'^; this 
value is comparable to what is observed in the Galactic 
HII regions or giant molecular clouds. To estimate the 
size of the neutral gaseous cloud, we adopt the observed 
log Af(Zn+) > 13.7 and assume a solar metallicity to infer 
log A^Hi > 21. We derive a cloud size of I = iVn / nn > 3 
pc. Lower metallicity would imply a higher A^hi and a 
larger size for the cloud. 

5. IMPLICATIONS FROM OBSERVATIONS OF 
FINE-STRUCTURE TRANSITIONS 

Additional constraints on the gas density, UV radiation 
intensity, and the spectral shape of the afterglow emission 
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can be learned from comparing the line strengths of dif- 
ferent fine structure lines (e.g. Bahcall & Wolf 1968). Al- 
though Fe II and Si II fine-structure transitions have yet 
to be observed in intervening quasar absorption line sys- 
tems (e.g. Howk, Wolfe & Prochaska 2005), they appear to 
be a generic feature in the GRB host environment. Their 
presence indicates unique physical conditions in the ISM 
hosting the GRBs. 

The excited states can be produced through either pho- 
ton pumping or coUisional excitations. In this section, 
we will investigate these various excitation mechanisms. 
Specifically, we will consider for photon pumping (1) di- 
rect excitation by IR radiation from sources local to the 
gas^ and (2) fluorescence following excitations by ultravio- 
let photons. For collisional excitations, we will discuss (3) 
collisions with charged particles, primarily electrons, and 
(4) collisions with neutral hydrogen. Although the photon 
pumping processes are generally weak in typical astrophys- 
ical environments, the radiation field of the GRB afterglow 
is extreme (Figure 3) and can dominate the excitation pro- 
cesses. 

To estimate the strength of the afterglow UV radiation 
field at the location of the absorbing gas in GRB 051111, 
for example, we adopt the minimum distance from § 4.1 
and Equation 3, r^-m — 80 pc. We conclude that the ma- 
jority of Fe"*" resides at a distance of 80 pc or greater from 
the afterglow. At r = 80 pc, the FUV radiation field is 
still > lO^G/Go (Figure 3), stronger than that observed 
for intense star- forming regions. 
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Average apparent optical depth profiles normalized by 



the degeneracy {gj = 2 J -|- 1) of the five levels of Fe+ for 

GRB 051111. One notes that the excited states track one another 
closely but that the ground state exhibits greater optical depth at 
D > km s"-*^ relative to the excited states. Also note that the gj- 
normalized column densities of the J = 9/2 ground-state level are 
scaled down by a factor of 20 to enable comparison with the excited 
states. 



Figure 5 shows the apparent optical depth profiles of 
the five a^D levels of Fe+ from GRB 051111 normalized 
by their degeneracy gj = 2J + 1. We emphasize three 

®Note that the CMB temperature and photon density at z = 
1.54948 cannot account for significant external IR pumping. 



observations from the Figure. First, the velocity profiles 
of the excited states (color histograms) track each other 
very closely. Second, the gj normalized profiles of the 
excited states have roughly the same value. Indeed, the 
total column densities of the J = 7/2toJ=l/2 states are 
roughly proportional to the degeneracy Nj oc gj with a 
modest but statistically significant decrease with increas- 
ing J. To populate the excited levels according to their 
degeneracy implies a saturated excitation process, i.e., an 
excitation rate that is much faster than the spontaneous 
decay rate. Therefore, it is unlikely the physical conditions 
vary considerably with velocity and we can integrate the 
profiles for the excited states over the entire velocity range 
to consider the total column densities. Lastly, the ground 
state {J — 9/2) exhibits a discrepant gj normalized profile 
from the excited states, containing « 20x higher column 
density than the value expected from the population of the 
excited levels. The ground level is clearly not populated 
according to its degeneracy relative to the excited states. 

Because no viable excitation mechanism will populate 
only the excited states according to gj, we conclude that 
the sightline has intersected at least two phases of gas: 
one with minimal excitation which contains roughly 95% 
of the ground-state gas and one where all of the a^D lev- 
els are roughly populated according to gj. These two 
phases could be embedded within one another (e.g. a dense 
core within a less dense cloud) or even represent sepa- 
rate layers of the same clump. Figure 5 offers additional 
kinematic evidence for two distinct phases; one observes 
that the ground-state profile is more sharply peaked at 
V = km s^^ than the excited states. For these reasons, 
in the following discussion we will not include observa- 
tional constraints drawn from the ground states of Si+, 
Fe+, or O". While Figure 5 presents only the results for 
GRB 051111, we note that the data from GRB 050730 
is qualitatively similar; the excited states are populated 
according to gj and the ground state shows a significant 
overabundance (Table 2). 

5.1. Direct IR Pumping 

We first examine the population ratios of the excited 
states under the scenario of excitation by IR photons. Di- 
rect IR pumping can populate the excited states of Si"*", 
Fe"*", and 0° through forbidden magnetic dipole transi- 
tions. Because these transitions have selection rules A J = 
0,±1, excited states beyond the first level must be popu- 
lated by an upward cascade of transitions. If IR pumping 
dominates, then the excitation of level j relative to the 
ground-state i is given by 



(6) 



N, ga + 

where is the number of photons at the appropriate 
frequency v 



N7. 



Sirhc 



(7) 



To explore the process of IR pumping, we consider the 
spectrum of the radiation field required to reproduce the 
observed column density ratios between the excited states 
of Fe+ for GRB 051111. We ignore the J = 9/2 state 
because « 95% of the ground-state ions must arise in a 
distinct region from the excited Fe^ gas. Figure 6 presents 
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Fig. 6. — Calculated intensity of the radiation field required to 
reproduce the relative populations of the J = 7/2,5/2,3/2 and 
J = 1/2 levels of Fe+ for GRB 051111 via direct IR pumping (mag- 
netic dipole transitions). The error bars reflect l-tr uncertainties 
in the ratios. For the upper two transitions (lower frequency), the 
observed ratios are within la of satisfying rugj = njQi. There- 
fore, there is no formal upper bound to ly and we show the central 
value but also express the result as a lower limit. The blue dashed 
curve qualitatively represents the flux from the GRB afterglow, i.e. a 
radiation field with oc v~^-^ . The red dotted curve shows the ra- 
diation field of a thermal source with In (x . The latter radiation 
field is ruled out by the observations. 



the intensities inferred from the observed population ratios 
of two adjacent levels (J = 7/2 -> 5/2, J = 5/2 ^ 3/2, 
and and J = 3/2 — > 1/2) versus the corresponding exci- 
tation frequency. The error bars reflect 1-a uncertainties 
in the ratios. For the upper two transitions, the observed 
ratios are within Icr of satisfying Uigj — n-jgi. Therefore, 
there is no formal upper bound to Ii, and we show the 
central value but also express the result as a lower limit. 

Superimposed on the plot are the predictions for two 
radiation fields with oc ly^ (e.g. a thermal spectrum 
with T > 200K) and oc (the optical afterglow 

spectrum; Butler et al. 2006). The thermal radiation field 
is ruled out at high confidence, whereas the GRB spec- 
trum (or a slightly steeper spectrum, /3 « 0, which might 
be appropriate for these low frequencies) is acceptable. 
For GRB 050730, where the population of the excited 
levels are all consistent with gj due to large measure- 
ment uncertainties, the data are consistent with a constant 



10 ^ ergs ^ 



cm 



Hz- 



over t/ = 2 - 9 X 10^^ Hz. 



This is higher than any reasonable IR source (i.e. dust) in 
even the most intense star-forming regions, and could only 
be produced by the GRB afterglow. 

Next, we estimate the number of Fe"'" ions that were ex- 
cited by IR photons from the afterglow prior to the start 
of the echelle observations. We calculate the excitation 
rate as a function of distance from the afterglow, taking 
into account the spontaneous decay rate of the excited lev- 
els. The spontaneous decay coefficient of the first excited 



level is A 



7/2,9/2 



— 2.12 X 10 s giving a decay time of 



8 minutes and an oscillator strength of 



Ajito 



2.16 X 10" 



(8) 



where to — 2.419 x 10~ s, a is the fine-structure con- 
stant, and Eji is the transition energy (Ryd). For a gas 
'cloud' with column density iV(Fe"'") = 10"'^^ cm~^ and 
Doppler parameter h — \b km s~ 
line center is 

1.5 X 10-2iVA/,-, 
To = ; 



the optical depth at 



5.6 X 10" 



(9) 



with all quantities expressed in cgs units. The estimated 
optical depth indicates that the gas is optically thin to the 
IR photons. To estimate the number of photons absorbed 
by the cloud per decay time, we calculate the equivalent 
width in the optically thin limit 

Wx = 8.85 (^^^JL_] ( -A_^ = 13mA . 



1015 cm- 



1000. 



(10) 

Figure 3 shows the number of IR photons emitted from 
tohs = 50 s to 3840 s with energy -£9/2,7/2 — 385cm~i in a 
lA window. The number of excitations per ion'' per decay 
time is 



JR 



lA 



iVion 



Aide 



At 



1 pc 



26 



pc 



= 10 



(11) 
ions 
1487 s. 



15 



where NJ — 10^^-^ photons per cm ^, Uio, 
per cm-2, Atjecay = 480 s, = Atobs/(l ■ 

According to Equation 11, there is a sufficient photon 
flux at r = Ipc from GRB 051111 to excite the Fell 
J = 7/2 level prior to the onset of the echelle observations 
but the photon flux is too low at r = 80 pc. Our anal- 
ysis therefore indicates that IR pumping is unlikely the 
primary mechanism for producing the excited ions found 
in GRB 051111. Furthermore, in the next sub-section we 
will show that the excitation rate by UV pumping is or- 
ders of magnitude higher than IR pumping. Therefore, we 
conclude that IR pumping is an unimportant excitation 
mechanism for gas observed in GRB afterglow spectra. 

5.2. Indirect UV Pumping 

The fine-structure states of Fe''", Si"*", and O'^ ions can 
also be populated by the absorption of a UV photon to an 
upper level, followed by spontaneous radiative decay to an 
excited level of the ground term. Requiring electric dipole 
transitions, the maximum angular momentum change is 
|AJma2;| = 2. Therefore, one requires a series of excita- 
tions to populate the J = 3/2,1/2 levels of Fe"*", and only 
a single absorption for the Si+ ion and O*^ atom. The anal- 
ysis of UV pumping is similar to IR pumping, although it 
requires a careful consideration of all the upper levels. To 
perform the calculation, we adopt the software packaged 
Popratio developed by Silva & Viegas (2001, 2002) which 
incorporates modern atomic data for these various ions. 

Similar to the treatment of IR pumping, we investigate 
(1) whether the observations describe a physically sensible 
radiation field; and (2) whether the pumping rate is suffi- 
cient to excite a significant column of gas prior to the start 
of spectroscopic observations assuming the GRB afterglow 
is the radiation source. We first examine the radiation 
field required to explain the observed population ratio of 
different excited states. 

'^Note that the rate is independent of the number of ions assumed 
here provided the optical depth remains small. 
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Fig. 7. — Indirect UV pumping of the first excited state relative to 
the ground state for Si+, Fe"*", and O" versus the far-UV intensity 
(normalized to Habing's constant, Go). At large intensity, the levels 
saturate to ratios according to their relative degeneracies gj. Note 
that Si"*" exhibits a much higher excitation ratio than the other gas. 



Figure 7 shows the excitation of the first excited states 
of Fe+ , Si+ , and 0° relative to the ground state as a func- 
tion of the far-UV intensity in terms of Habing's constant 
Go = 1-6 X 10"'^ erg cni~^ s~^. The resuhs are drawn 
from Popratio calculations assuming physical character- 
istics for the gas such that UV pumping dominates (i.e. 

— 1 cm~'^). We emphasize three points from the figure. 
First, all of the excited states saturate at G/Gq > 10^. At 
these intensities, the levels are populated only according to 
the degeneracy of the excited state relative to the ground 
state and the differences between ions reflect the various 
J values. Second, the excitation fraction is negligible and 
is nearly proportional to G for G/Gq < 10^, beyond which 
pumping induced excited ions become dominant. Finally, 
Si"*" exhibits « 10 times higher excitation fraction than 
Fe"*" and for the same G/Gq values. This is partly 
due to the fact that the Si II levels are inverted {gj for 
the excited state is twice that of the ground state), but 
also arises from the detailed atomic physics of the upper 
levels^. 

If we adopt a minimal excitation for Fe+ , then for 
GRB 051111 only roughly 7% of the observed 7Vj=9/2(Fc+) 
originates in the excited phase, then our observations im- 
ply a UV radiation field of G/Gq > 10^ under a steady- 
state assumption (Figure 7). At this radiation intensity, 
we also infer a population ratio for Si"*", ?^J=3/2/'^^ .7=1/2 ^ 
2, indicating that roughly 5% of the observed ground-state 
Si'^ originates in the excited phase. This exercise shows 
that roughly the same fraction of the observed Fe+ and 
Si+ ions have been excited toward GRB 051111 under the 
indirect UV pumping mechanism. 

Examining the far-UV radiation field of a GRB after- 
glow (Figure 3), we note that the unattenuated intensity is 
sufficient in steady-state conditions to excite all of the ions 
out at even several hundred parsecs where G/Gq > 10^. 



"It is not an effect of UV absorption strengths; all of the ions 
exhibit a significant number of resonance transitions in the far-UV. 



In particular, the luminosity of the GRB exceeds that for 
star-forming regions with SFR=1OOM0 yr~^ by a factor 
of 1000. Therefore, if the excitation rate of the ions is 
large on the time-scale of the afterglow observations (e.g. 
minutes to hours), then the burst should excite Fe"*", Si"*", 
and 0° via indirect UV pumping. 

The excitation rate due to fluorescence transitions can 
be roughly estimated under the assumption that the gas 
is optically opaque to many of the resonance transitions. 
In contrast to IR pumping where the process occurs in 
a cascade of magnetic dipole transitions, UV pumping 
proceeds through many transitions (~ 10 — 20) from the 
ground-state multiplet to upper- level multiplets. A de- 
tailed radiative transfer calculation is necessary given the 
transient nature of the GRB. Such a calculation is beyond 
the scope of this paper and will be the subject of a future 
work (Mathews & Prochaska, in prep). Here we present a 
simple framework for such calculations. 

To estimate the number of far-UV photons absorbed 
per decay time, we first consider Si"^. Adopting a col- 
umn density iV(Si~*') = lO^^'^cm"^ and Doppler param- 
eter b = lOkms^^ which roughly matches the profile 
of the excited gas for GRB 051111, we first produce a 
model spectrum of the ground-state Si II transitions with 
912 A < Arost < 2000 A assuming the atomic data given 
in Morton (2003). We then measure the equivalent width 
of each si II transition which allows us to calculate the 
number of photons absorbed per second from the GRB af- 
terglow spectrum (Equation 1). Summing over all of the 
transitions and evaluating the GRB afterglow at tobs = 50s 
and r — 10 pc, we calculate the total number of photons 
absorbed per minute per ion 

-2 / ,r \ -1 



80 



5.4 X 1031 



10 pc 



N,, 



10 



14.5 



cm 



assuming the unattenuated spectrum given by Equation 1. 
Although only a fraction of these photons will excite the 
Si+ ions (roughly 66%), the rate is sufficiently high to show 
a modest Si II* column density at r = 50 pc at < = 20 min 
after the burst. Even at r = 100 pc, there would be w 16 
excitations per ion prior to the onset of observations for 
GRB 051111. Therefore, we consider the excitation rate to 
be sufficient out to at least r = 100 pc from this afterglow. 

We carry out similar calculations for O" and Fe+ and 
find results that match Equation 12 to within a factor of 3 
(Fe"*" has a higher rate and 0° is lower, primarily reflecting 
the rest wavelengths of the resonance transitions). We em- 
phasize that these excitation rates exceed IR pumping by 
more than an order of magnitude (compare Equations 11 
and 12). Therefore, we conclude that indirect UV pump- 
ing is the principal excitation mechanism in GRB when 
photon pumping dominates. 

We note that the excitation rate per ion given by Equa- 
tion 12 is inversely proportional to the number of ions 
assumed. This is because the gas is optically thick in the 
majority of Si II transitions. The same holds for Fe II and 
O I transitions. Therefore, the equivalent width (hence, 
the number of photons absorbed) is insensitive to the col- 
umn density assumed provided Nion ^ lO^^cm"^. Given 
the transient nature of the GRB, this has important im- 
plications for the excitation of these ions. Most impor- 
tantly, for a gas cloud (or clouds) with total column den- 
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sity Nion > 10^^ cm"^, only a fraction of the gas along the 
sightline is excited prior to the onset of observations. 

To zeroth order, the number of excited ions is given by 
the column density where the majority of transitions reach 
the saturated portion of the curve-of-growth. We have in- 
vestigated this column density for the Fe"*", Si"*", and O*' 
gas. We find that for 6 = 10 km s^^, the excitation sat- 
urates at A'^f™ « 10^*-^ cm~^ for Si+ and Fe+ and a few 
times larger value for 0°. To first order, this limiting 
column density scales with the Doppler parameter of the 
gas, i.e. wider profiles should show larger column densi- 
ties of excited gas when the excitation rate is limited by 
line saturation. We suspect that is not a coincidence that 
the column densities of excited gas in GRB 051111 and 
GRB 050730 are roughly consistent with this calculation. 

In summary, UV pumping is a viable mechanism for 
excitations of the fine-structure states observed in GRB 
afterglow spectroscopy. For the case of GRB 051111, wc 
find the gas would be pumped to at least 100 pc from the 
GRB afterglow if the far-UV light is not heavily attenuated 
by dust or intervening metal-rich gas. We draw similar 
conclusions for GRB 050730 and the other GRB listed in 
Table 3. In short, the radiation field of GRB afterglows 
is likely to excite 81+ to at least 100 pc in the majority of 
systems observed to date. 



1.0 




Fig. 8. — Indirect UV pumping of the Fe"*" excited levels relative 
to the J = 7/2 level (normalized by the degeneracy of the states gj) 
as a function of the far-UV intensity. The levels saturate at high 
intensity, i.e., they are populated according to gj. Also note that 
the J = 3/2 and J = 1/2 levels track one another at all intensity. 
This occurs because the two states require the absorption of two 
photons to be populated via UV pumping. The signs indicate the 
observed ratios of the excited levels for the gas toward GRB 051111. 
These observations are consistent with an intensity G/Gq ~ 10^'^. 



from the trivial case where the levels are saturated corre- 
sponding to Tex ^ Ej^i/2/k. Indeed, this is the case for 
GRB 050730. In Figure 8, we present the predicted abun- 
dances of the J = 5/2,3/2 and 1/2 levels relative to the 
J = 7/2 state normalized by the degeneracies gj as a func- 
tion of the far-UV intensity. In contrast with IR pumping, 
the spectral shape is unimportant to this calculation; the 
results only depend on the overall excitation. 

An important prediction for UV pumping is that the 
J = 3/2 tracks the J — 1/2 level for all intensities. At 
G/Gq > 10^, the levels are saturated, i.e. populated ac- 
cording to gj. At lower intensity G/Gq < 10^, the relative 
populations are strictly inconsistent with a Boltzmann dis- 
tribution. For intermediate values {G/Gq ~ 10^'^), how- 
ever, the differences from a Boltzmann distribution are 
small. The plus signs marked in Figure 8 show the ob- 
served ratios for GRB 051111 and indicate the data are 
well fit by a pumping scenario with a single intensity. 

To conclude, we find that UV pumping can reproduce a 
Boltzmann distribution of Fe II levels with Tex > 2000K. 
A key prediction of UV pumping is that the relative pop- 
ulation of the J = 3/2 and 1/2 levels is given by the ratio 
of their degeneracies: '^3/2/'^i/2 — 2. This is inconsistent 
with a Boltzmann distribution for Tex < 1000 K and is 
a direct test for distinguishing between coUisional excita- 
tion and UV pumping: the observation of 113/2 /ni/2 > 2 
would rule out UV pumping. In the cases of GRB 051111 
and GRB 050730, however, we find that n3/2/'T-i/2 ~ 2 
and that both UV pumping and coUisional excitation are 
viable mechanisms. 
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Fig. 9. — CoUisional excitation of the first excited states relative 
to the ground state for Si+, Fe+, and versus the electron density. 
Note that for O", neutral collisions dominate and we have assumed 
an ionization fraction x = 10~^. 



While we have demonstrated that GRB afterglows are 
sufficiently bright to UV pump gas to large distance on 
short time-scales, we must also consider whether the rel- 
ative populations of the excited states reproduce the ob- 
served ratios. In the next section, we will demonstrate that 
the Fe II excited levels for GRB 051111 are well described 
by a Boltzmann distribution with excitation temperature 
Tex = 2600K. Because UV pumping is not a thermal pro- 
cess, one does not expect a Boltzmann distribution aside 



5.3. CoUisional Excitation 

The final excitation mechanism to consider is coUisional 
excitation. For the Fe"*" and Si"^ ions, collisions with elec- 
trons are likely to dominate the rates^. For 0°, the rates 

^We caution that the neutral collision rates for Fe+ have not yet 
been precisely calculated or measured in the laboratory. 
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for collisions with hydrogen atoms are much faster and are 
dominant for a mostly neutral cloud, i.e., x < 10^^. The 
collision rates do vary with temperature but by less than 
an order of magnitude for temperatures T = 100 to 10000 
K. In the following, we will present results for T — 2600 
K. 

In Figure 9 we present the population the first excited 
state relative to the ground state as a function of electron 
density Ue (neutral hydrogen density ?ihi for O''). The 
results are qualitatively similar to those for UV pumping 
(Figure 7) where one exchanges density with far-UV in- 
tensity. At large density, the excited states are populated 
according to their relative degeneracies although reduced 
by the Boltzmann function (see Equation 13 below). At 
low density, the population is proportional to the density 
of the colliding particle. Also, the excitation of Si"*" oc- 
curs at significantly lower density although note that the 
excitation of O'' would be comparable for gas with low 
ionization fraction. In contrast to UV pumping, the ex- 
citation of Fe"*" requires significantly higher densities than 
Si"'" and O'^. One finds Ue > 10^ cm^'^ for even a modest 
population of the J = 7/2 level. 

When the collisional de-excitation rate qjiUe exceeds the 
spontaneous decay rate Aji, the excited states are popu- 
lated according to a Boltzmann distribution, 

!!i= 5^exp[-(ii;,,/fc)/TE,] , (13) 
rij gj 

where Eij = Ei—Ej is the difference in energy between the 
two states, and Tex is the excitation temperature. Com- 
paring predictions following Equation 13 with the observed 
column densities of ions at different excited states allows 
us to examine whether collisional equilibrium is a represen- 
tative model for the excitation mechanism. A consistent 
fit in turn leads to a strong constraint for the gas temper- 
ature. 




400 



600 



800 1000 
E,/k (K) 



1200 



1400 



Fig. 10. — Plot of Ni/gi for the Fe+ excited states toward 
GRB 051111 as a function of their energy Ei above the ground state. 
Overplotted on the figure is the function AeK'p[—{Eij /k) /Tex for a 
minimum value of T^x = 2700 K. The inset shows the reduced 
values for the best fit for a range of T^x values. 



Figure 10 shows for GRB 051111 the observed column 
densities scaled by the corresponding degeneracy Ni/gi for 



the Fe II excited states as a function of the energy Eij 
above the ground state J = 9/2. The error bars reflect 
l-CT uncertainty in Ni . The modest decrease of Ni j gi with 
Eij suggests the levels are Boltzmann populated with an 
excitation energy greater than the largest Eij value. The 
solid (blue) curve indicates the best-fit Boltzmann func- 
tion with a best-fit excitation temperature Tex = 2600 K. 
The inset shows the minimum ^ values as a function of 
Tex; the 99% c.l. limits are rEx(99%) = 1800 - 6000 K. 
The reduced ^ is nearly unity, = 0.99, supporting the 
assumption that the Boltzmann function is a representa- 
tive model. 
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Fig. 11.— Same as Figure 10 except for GRB 050730. 



We have repeated this analysis for GRB 050730 and the 
results are presented in Figure 11. The column density 
measurements have larger uncertainty and therefore place 
weaker constraints on Tex- Nevertheless, the Fe+ excited 
state population is well fit by a Boltzmann function with 
Tex ~ 6000K with a firm lower hmit of lOOOK. 

5.4. Summary 

After examining the various mechanisms for exciting 
Fe II, Si II and O I excited levels in the presence of a 
GRB afterglow, we have demonstrated that indirect UV 
pumping and collisional excitation are both viable mech- 
anisms for producing the excited ions. For indirect UV 
pumping to work, it requires an intensified far-UV radia- 
tion field: G/Gq > 10^. For collisional excitation to dom- 
inate, it requires a dense environment of > 10"^cm~^. 
We rule out IR pumping for gas at distances > 10 pc and 
also find that the excitation rate is small compared to UV 
pumping at all distance (compare Equations 12 and 11). 
Therefore, the excitation of gas near GRB afterglows is 
primarily determined by indirect UV pumping and/or col- 
lisional excitation. 

Another important conclusion is that the predictions for 
UV pumping and collisional excitation are very similar. 
For example, there is the trivial case at very high intensity 
or density where the excited levels are principally popu- 
lated according to their degeneracy gj, modified by the 
Boltzmann factor collisions. Although the detailed level 
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Table 4 

TESTS OF UV PUMPING AND COLLISIONAL EXCITATION 



ObbC'ivaLioii 



JV^=3/2(Fe+) = 2iV^=i/2(Fe+ 
Inverted Fe+ population 
Variability 

Absence of atomic gas 

Overabundance of the Fc+ J - 
[A^J=7/2(Fc+)-W;=3/2(Si+) 

Arj=i(O0)<7Vj=3/2(Si+) 

Nexcited > 10^^ cm 2 



IJ\' i filliping 



(.'oliisiuiis 



= 9/2 level 
< [-1 + e(Fe/Si)]'' 



Required 
Ruled out 
Expected 
Expected 
Expected 
Ruled out 
Ruled out 
Pre-burst only 



rie > lO** cm 
Ue < 10* cm-3 
Ruled out 
Challenging 
Challenging 
Allowed 

Allowed for X > 10 
Allowed 



T > 2500K 



"e(Fe/Si) is the intrinsic, logarithmic Fe/Si gas-phase abundance allowing for nucleosynthetic and differential depletion effects. 



populations arc different because of the Boltzniann factor, 
one requires very high SNR observations to distinguish the 
cases if Tex > 2000K. In the next section, we will consider 
a few additional tests for distinguishing between the exci- 
tation mechanisms. 

Regarding the gas near GRB afterglows, it is clear from 
Table 3 that photoionization dominates at distances less 
than 10 pc from the GRB and radiative processes are sig- 
nificant to 100 pc or further for the majority of afterglows. 
We conchidc that circumstcUar gas associated with the 
GRB progenitor will be excited by UV photons, if the gas 
is not entirely ionized by the event. Therefore, it is dif- 
ficult to determine the density and temperature of such 
gas from the analysis of fine-structure transitions. A full 
calculation of the optical depth and time dependent pro- 
cesses of excitation are warranted and could help reveal the 
physical conditions in the absorbing medium (Mathews & 
Prochaska, in prep). 

6. DISTINGUISHING BETWEEN PHOTON PUMPING AND 
COLLISIONAL EXCITATION 

In the previous section, we argued that the excitation 
of fine-structure states of Fc+, Si+, and gas can occur 
along the sightline to GRB afterglows because of indirect 
UV photon pumping and/or coUisional excitation. Each 
process has different implications for the physical condi- 
tions of the gas (e.g. density, temperature). It is therefore 
important to distinguish between the two mechanisms. We 
summarize a set of tests in Table 4. 

There is a sufficiently large number of Fe II transitions 
covering a wide range of UV wavelengths in the rest frame 
and nearly every GRB afterglow spectrum includes the full 
set of Fe II excited levels (c.f. GRB 050730 and GRB 051111) 
Those sightlines which exhibit positive detections of the 
excited states allow different tests of the excitation mech- 
anism; we consider the two cases separately. 

6.1. When Excited States of Fe'^ are Detected 

We have shown that under excitation equilibrium the 
indirect UV pumping scenario predicts a population ratio 
that is proportional to the degeneracy of the excited level, 
while the collisional excitation scenario predicts a popula- 
tion ratio that follows the Boltzmann function. In prin- 
ciple, one can therefore distinguish between these two ex- 
citation mechanisms based on the observed relative abun- 
dances between different excited ions. Our analysis has 



shown, however, that measurement uncertainties often pro- 
hibit us from achieving the goal. Despite the caveat in 
obtaining more precise measurements, we highlight in this 
section two features that supports photon pumping and 
two features that supports collisional excitation for future 
studies. 

Under photon pumping, we first note in Figure 8 that 

the J = 5/2 state is significantly more populated than the 
higher levels at low intensity G/Gq < 10^. This can be 
understood by the fact that one can populate the J = 5/2 
and J = 7/2 states with a single photon whereas excita- 
tion to the J = 3/2 and J = 1/2 levels requires at least 
two-photon processes. Second, detection of line variabil- 
ity offers another discriminator between UV pumping and 
collisional excitation processes. For collisional excitation, 
we do not expect variability on time-scales that could con- 
ceivably be probed by GRB afterglows (i.e. less than a 
few days). In contrast, the absorption- line strength is ex- 
pected to vary under photon pumping, because (1) Fe^ 
becomes more ionized with time due to the afterglow; (2) 
the afterglow radiation field declines with time, leading to 
decreasing G/Gq; and (3) the diffusion of photons through 
an optically thick cloud is a time dependent process and 
could lead to both higher and lower column densities of the 
excited levels (Mathews & Prochaska 2006). To observe 
both features require that additional echelle observations 
be carried out at a later time. 

We also note that in one extreme there is a maximum 
column density N!^^^ of ions that can be excited by UV 
pumping prior to the onset of spectroscopic observations, 
because of optical depth effects. The value of N^j^^ is pri- 
marily established by the velocity structure of the medium 
(i.e. the relative velocity of the absorbing clouds and their 
Doppler parameters) and to a lesser extent the size of the 
cloud (1/r^ dimming). An accurate estimate of A^moL will 
require a time dependent radiative transfer calculation, 
but we note that nearly all of the Fe II, O I, and Si II 
transitions saturate at N^^^ « lO^^cm"^. Therefore, if 
we observe a column density of excited gas significantly 
in excess of lO^^cm"^, this would be strong evidence for 
collisional excitation, or possibly pre-burst UV pumping 
as optical-depth effects would not prevent the entire cloud 
from being excited (Sarazin, Rybicki & Flannery 1979). 

Under collisional excitation, we first note that at the im- 
plied gas density rie > 10^ cm~^ (Figure 9) we also expect 
to observe large fractions of the elements C, Fe, and Mg 
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in their atomic states. Although we have shown in § 4.1 
that the GRB afterglow can photoionize atomic gas at 
large distances, there is a relatively narrow window where 
the gas would be ionized and coUisional excitation would 
dominate indirect UV pumping. Therefore, the presence 
of Fe°, Mg°, or C" may be viewed to support coUisional 
excitation. 
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Fig. 12. — CoUisional excitation of the Fe II excited levels relative 
to the J = 7/2 level normalized by the degeneracy of the states 
Pj as a function of the electron density. The levels saturate at 
high density, i.e., they become populated according to gj weighted 
by the Boltzmann factor exp{—AE /kTEx) where we have assumed 
Tex = 2600K. Also note that at intermediate density the relative 
populations are roughly the same and that at low density the J = 
1/2 and J = 3/2 exhibit an inverted population. 



Second, we note in Figure 12 that under coUisional exci- 
tation for a gas temperature Tex — 2600 K the population 
ratios between different excited levels are significantly dif- 
ferent from the Boltzmann distribution at < lO'^cm"'^. 
Specifically, both J = 5/2 and 3/2 levels are significantly 
underabundant relative to the J — 7/2 level. In addition, 
the population becomes inverted, especially the J — 1/2 
level. The relative abundance at these lower densities con- 
trast with anything predicted for UV pumping. Although 
Figure 12 applies to a single phase of gas with T = 2600 
K, the results are similar for any temperature T < 10000 
K. It is, however, very challenging to observe an inverted 
population of Fe II levels because the absolute excitation 
rate of Fe+ is small for rie < lO'^cm^'^. 

An intriguing aspect of the Fe II excited levels for GRB 0511 
and GRB 050730 is the overabundance of the ions in the 
ground level. We also observed for GRB 051111 that the 
ground-state of Si"*" is at least 10 x more populated than 
the excited state. While one could construct a scenario 
which reproduces the observations for either excitation 
mechanism, we interpret the overabundance of the ground 
level as evidence for UV pumping for two reasons. First, 
the effect is naturally explained within the context of UV 
pumping by optical depth effects. If the cloud is optically 
thick, then one expects partial excitation especially if one 
cloud 'shadows' the rest of the gas that is further away 
from the afterglow. Second, the similarity in kinematics 



suggest the gas is co-spatial and yet under coUisional ex- 
citation the excited gas would require far greater pressure 
than is conceivable for the unexcited phase. 
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Fig. 13. — Comparison of the predicted column densities of the 
first excited states of Fe+ (black) and (blue) relative to the 
J = 3/2 excited level of Si+ for indirect UV pumping (solid curves) 
and coUisional excitation (dotted; T = 2600K). We have assumed 
gas-phase abundances e(0/Si) = 1 and e(Fe/Si) = 0. It is important 
to note that the coUisional excitation calculations for O'^ assume an 
ionization fraction x = 10~^. Higher values of x would lead to lower 
N{Ol*) values relative to A''(SilI*) in a roughly linear fashion. 



6.2. When Excited States of Fe^ are not Detected 

For GRB afterglows where the excited Fe II levels are 
not detected, the abundances of the first excited levels of 
O" and Si"*", together with an upper limit in excited Fe"*" 
can be applied to constrain the excitation mechanism. Fig- 
ure 13 presents the values of Nj=i{0^) and A^,7=7/2(Fe"'") 
relative to iV,/=3/2(Si"'") as a function of electron density 
and intensity for coUisional excitation (dotted curves) and 
UV pumping (solid curves), respectively. Note that the 
curves assume that the gas-phase O" abundance is 10 times 
higher than Si+, e(0/Si) = logAf(0°) - logiV(Si-^) = 1, 
and that Si"*" and Fe+ have equal gas-phase abundance, 
e(Fe/Si) = 0. Differential depletion, in particular, would 
modify these assumptions: increasing e(0/Si) and decreas- 
ing e(Fe/Si). 

We first focus on Fe+ and Si^. It is evident that col- 
lisional excitation at n,, < 10^ allows for significantly lower 
lvalues of A^j^7/2(Fe+)/A^j^3/2(Si^) than UV pumping with 
G/Gq < 10^. The main result is that UV pumping does 
not allow log 7Vj^7/2(Fe+) -log Arj^3/2(Si+) < [-l+e(Fe/Si)]. 
Therefore, a GRB afterglow spectrum exhibiting Si II* ab- 
sorption and strict upper limits on the excited Fe II tran- 
sitions is likely to indicate coUisional excitation. There 
are, however, two caveats. First, it assumes knowledge of 
the inherent Fe/Si ratio. In general, this may be inferred 
from the ground state populations. Second, the results in 
Figure 13 also assume steady-state conditions. The obser- 
vations of GRB afterglows may occur on time-scales that 
are short compared to the lifetime rufo of the excited lev- 
els (especially for Si+ and O" which have rufe = 1.3 and 
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3.1 hr, respectively). Therefore, a significant underabun- 
dance of A^j=7/2(Fe+) relative to iVj^3/2(Si"'") is possible 
from UV pumping if one observes the excited gas at a time 
when the excited Fe+ ions (riifc « 8min) have decayed to 
the groimd-state. In this case, of course, one would ex- 
pect line-strength variability between early- and late-time 
spectroscopic data. 

Next, we consider the and Si"^ predictions. Indepen- 
dent of the excitation mechanism, we find Nj^i{0'^) > 
-^j=3/2(Si''') for the assumed ionization fraction-^", x w 
rie/nn = 10~^. For a predominantly neutral gas, coUi- 
sional excitation of 81+ will roughly match that of O'' and 
^./=i(0") > A^j=3/2(Si+) is predicted. For larger values 
of X, however, the excitation rate for 81+ is considerably 
higher and a significantly smaller ratio 
of A^,7=i(0")/iV,/=3/2(Si+) is expected. In fact, the only 
process which allows iVj=i(0°) < A^j=3/2(Si"*') is coUi- 
sional excitation in a predominantly ionized medium, e.g. 
the Carina nebula in the Galaxy (Walborn et al. 2002). 

7. DISCUSSION 

7.1. Implications for GRB 051111 and GRB 050730 

Despite different spectral coverage in the rest frame of 
GRB 051111 and GRB 050730, the echelle spectra of these 

afterglows exhibit ground-state 81"*" and Fe"*" features from 
all of their excited levels. In the case of GRB 050730, ex- 
cited O" is also present. But due to measurement un- 
certainties, their population ratios may be explained by 
Boltzmann distribution with Tex = 2600 - 6000K (Fig- 
ures 10 & 11), as well as by UV pumping provided G/Gq > 
10^'^. We have investigated variations in the column den- 
sities and set upper limits to the line-strength variability 
of < 10% per 1800s for GRB 051111 and < 30% per 1800s 
for GRB 050730. It appears difficult to unambiguously 
rule out either UV pumping or coUisional excitation as 
the dominant excitation mechanism. The only strong ar- 
gument against coUisional excitation is the absence of Fe" 
and C° atomic gas for GRB 051111 and GRB 050730, re- 
spectively. 

Under coUisional excitation, the excited ions occur ir- 
respective of the presence of the afterglow and the gas is 
required to have rie > 10^ cm'^ and T ~ 2000 - 6000 
K to populate the Fc II excited levels. If the excited 
gas arises from a layer of high density (rie > 10^ cm""^), 
then based on the absence of Fe I we can place an upper 
limit on the ratio of the electron density to far-UV inten- 
sity He X (G/Go)"^ (Equation 5). Adopting the best-fit 
Boltzmann distribution, we calculate the total Fc+ col- 
umn density NtotFe~^ = 10^'*'^cm~^, including all levels. 
Adopting a conservative upper limit to the atomic state 
7V(Fe°) < 10^2 cm-2 (Table 1) which accounts for con- 
tinuum and statistical uncertainty, we therefore find from 
Equation 5 that x (G/Go)"^ < 1 for T < lOOOOK. For 
He > 10^cm~3, we find that G/Gq > 10^. Therefore, we 
conclude that we would only observe significant popula- 
tions of the excited Fe II levels for gas in the presence of 
a strong far-UV radiation field. 

To satisfy G/Gq > 10^, it requires the gas to lie within 
« O.lpc of a T = 40000K O star with R = lOi?©. While 

^'^This is a rough estimate of the ionization fraction for a predom- 
inantly neutral gas. 



the GRB progenitor star is expected to reside in a star- 
forming region and therefore near other massive stars, we 
consider it highly unrealistic that the gas is at > 100 pc 
from the GRB but also within 0.1 pc of an O star. 

A similar conclusion can be drawn for GRB 050730 using 
C I. We infer a C"*" column density by scaling A^(Fe-l-) = 
lO^^'^cm"^ from the excited region with solar relative 
abundances: A''(C"'") = 10^^-^cm~^. The observed upper 
limit N{C°)/N{C+) < 10-2-3 requires G/Go > 5 x lO"*, 
again assuming a conservative limit to the temperature 
(T < lOOOOK). For GRB 050730, we have no significant 
constraint on the distance of the gas from the afterglow. 
But, if coUisional excitation dominates indirect UV pump- 
ing (and that < lO'^cm^'^), then this gas must also be 
experiencing a strong local source of far-UV radiation. We 
conclude based on the null detections of atomic gas that 
indirect UV pum.ping by GRB afterglow radiation is the 
dominant excitation process. 

We also note that there are additional challenges to a 
coUisional excitation scenario. First, the gas would require 
a very large heating source to maintain the inferred tem- 
peratures for coUisional excitation. For the implied physi- 
cal conditions (e.g. T = 2600 K and > 10^ cm-^), [Fe II] 
cooling dominates. Using the coefficients given by HoUen- 
bach & McKee (1989), we calculate a rate of nAjpeii] = 
4 x IQ-^^ergs s"-'^ H"^ assuming 1/10 solar metallicity and 
that 90% of the iron is locked into dust grains. This implies 
a coofing rate tcooi ^ kT / {uA^-p^i^) < 500yr~-^. Including 
the cooling by meta-stable states of O I would imply an 
even shorter cooling time. Second, gas with T « 3000 
K is generally not stable to thermal perturbation. In the 
case considered here, the [Fe II] cooling curve is relatively 
insensitive to temperature variations but very sensitive to 
density changes. For example, an increase in nn increases 
n A [Foil] and drives the gas to higher density and lower 
temperature imtil T <C 1000 K. Third, the implied pres- 
sure nT is extremely high. One would require a strong 
mechanical process (e.g. shock) to confine the gas and it 
would be surprising to therefore find the quiescent kine- 
matic characteristics that are observed (Figures 1,2). 

Adopting UV pumping as the the excitation mechanism, 
we can constrain the distance of the gas from the GRB 
based on known afterglow radiation field. For GRB 051111, 
the gas must arise at r w 200 pc. This distance is large 
enough to prevent the ionization of Mg° but also close 
enough that the Fe II and Si II levels would be UV pumped 
and also show the relative populations that are observed. 
In the case of GRB 050730, where the excited Fe II levels 
are populated according to their degeneracies gj, the gas 
must lie within « 100 pc of the afterglow (Table 3). The 
ionizing column of photons at Ipc is 6 x lO^^cm^^. The 
gas is therefore not likely at < 3 pc, where the far-UV ra- 
diation field from the GRB afterglow will photoionize all 
of the C+, Si+, 0°, and Fe+ gas. 

7.2. Implications of Fine-Structure Lines in Other GRB 

Afterglows 

The detection of 8i II* transitions is nearly ubiquitous 
in GRB afterglow spectra (Vreeswijk et al. 2004; Chen 
et al. 2005; Berger et al. 2005). This stands in stark 
contrast with intervening quasar absorption line systems 
where only C I* and C II* fine-structure lines have been 
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observed (e.g. Srianand, Petitjean & Ledoux 2000; Wolfe, 
Prochaska & Gawiser 2003; Howk, Wolfe & Prochaska 
2005). To date, the observation of Si II* absorption for 
GRB have been interpreted as the result of collisional ex- 
citation and therefore evidence for large volume densities 
in the gas. We have demonstrated in this paper, however, 
that UV pumping is the most Ukely excitation mechanism 
if the gas is located within a distance of « 100 pc. This 
conclusion holds for all of the GRBs with reported Si II* 
detections (Table 3) except for GRB 050408 (Foley et al. 
2005). And, for sightlines exhibiting Fe II* absorption, we 
contend UV pumping is likely the only mechanism. 

In summary, we conclude that in the presence of an 
intensified UV radiation field from the afterglow indirect 
UV pumping alone offers a simple explanation for all the 
excited fine- structure transitions observed in GRB host en- 
vironment to date without invoking extreme ISM proper- 
ties, such as high gas density. If we take the gas den- 
sity allowed by the observed relative abundances of atomic 
species, uh = 100 — 1000 cm~^ (§ 4.2), the observed large 
^Hi (e.g. Chen et al. 2005) implies a cloud thickness of 
3 — 30 pc which is typical of giant molecular clouds (Turner 
1988). In lieu of diagnostics from the fine-structure popu- 
lations, one must rely on comparisons of the atomic states 
with higher ionization levels to constrain the gas density. 
We contend that the principal quantity inferred from the 
fine-structure lines is the distance to the GRB afterglow. 
At present, there is no obvious diagnostic for the gas tem- 
perature. 

7.3. Implications for Gas Within lOOpc of GRB 
Afterglows 

Our analysis reveals a number of implications for the 
study of gas located within 100 pc of the GRB afterglow. 
The first is that this gas may be identified by the absence 
of atomic absorption, in particular Mg I gas. Although 
the sightline is likely to penetrate Mg I gas at distances far 
from the GRB (e.g. within the galactic halo), differences in 
the linc-of-sight velocities should distinguish these clouds. 
Therefore, a study of the kinematic differences between 
Mg I and Mg II profiles in GRB spectra is warranted. 

In the absence of atomic gas (e.g. Mg I absorption) 
that coincides kinematically with the low-ion species, the 
gas may have small distance from the GRB afterglow and 
therefore may be circumstellar material from the progeni- 
tor. Unfortunately, the afterglow radiation washes out the 
most readily available diagnostics for studying the density 
and temperature of this gas. We note, however, that self- 
shielding of UV photons matching various excitation ener- 
gies may be significant, especially for the majority of other 
GRB sightlines where the fine-structure column density is 
signficantly lower than the ground-state (e.g. Vreeswijk et 
al. 2004; Savaglio & Fall 2004). We expect this may also 
be the case for GRB 051111 and GRB 050730. It is there- 
fore possible that a significant portion of the gas is not 
excited by the afterglow and we could infer physical con- 
ditions for the bulk of the medium. 

Another implication is that one must consider the fine- 
structure levels in the chemical abundance analysis. This 
will be especially important for sightlines with low total 
metal- line column density. The correction to A^(Fe"'") in 
GRB 050730, for example, is 33% and we estimate the 



correction to -/V(Si+) is greater than 50%. If one cannot 
account for the excited-states in the analysis (e.g. because 
the lines are too saturated/blended for accurate measure- 
ment), then analysis of S+ is preferred because it does not 
have fine-structure levels near the ground-state. 

Finally, we emphasize again that to observe gas at 10 pc 
one requires a large pre-existing column density of neutral 
or partially ionized gas prior to the GRB event. This in- 
cludes highly ionized species like C IV. Consider, as an ex- 
ample, the relatively spectacular C IV gas observed toward 
GRB 021004 which several authors have interpreted to 
arise from gas associated with the GRB progenitor (Schae- 
fer et al. 2003; Mirabal et al. 2003; van Marie, Langer 
& Garcia-Segura 2005; Fiore et al. 2005). At tabs = 1 
hr after the GRB, we estimate an H I ionizing column 
density of 3 x 10^^ photons cm~^ at 10 pc assuming the 
afterglow properties given in Table 3. This column sig- 
nificantly exceeds the A^hi value observed for the gas to- 
ward GRB 021004, therefore we assume that all of the 
low-ion gas within 10 pc was either ionized prior to the 
GRB or prior to the observations. This helps explain the 
absence of strong low-ion absorption^ ^. More importantly, 
the column density of C^^ ionizing photons with energy 
hv = 64.5 to 70eV at lOpc is f« lO^Ocm'^. If this flux 
is unattenuated by gas within 10 pc then it would signif- 
icantly ionize all of the C+^ ions within a few tens of pc 
from the GRB. At these energies, the principal source of 
opacity is IIe+, but we estimate that this gas will also 
have been burned away by the GRB and its progenitor. 
We contend, therefore, that the majority of C+^ gas ob- 
served along the GRB 021004 sightline is either associated 
with the larger star forming region and/or halo gas from 
the GRB host galaxy and/or its galactic neighbors. 
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APPENDIX 

PHOTOIONIZATION EQUILIBRIUM 

Under steady-state conditions and photoionization equi- 
librium, the balance of photoionization and recombination 
for Mg° and Mg+ gives 

^^The observation of Al II absorption is an interesting exception. 
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n{Mg°)n^ca^^^ = n(Mg+)ne [a.(T) + adi(T) ] (Al) 

where and a^i are the radiative and dielectronic recom- 
bination coefficients and a^ff is the cross-section to pho- 

toionization integrated over the incident radiation field. 
For the recombination rates, we wiU use the fitting for- 
mula given by Aldrovandi & Pequignot (1973), 

afs+(T) = 1.4 X 10-i3(y/i04)-°-855cm3s-i , (A2) 
and ShuU & van Steenberg (1982), 



a^^*{T) = 4.49 X IQ-^T-^-^ exp(-5.01 x 10'^ /T) [l +0.0021 exp(-2.61 x lOV?") ] cm^ s"^ 

For the photoionization rate, we integrate the Galactic far- 
UV radiation field determined by Gondhalekar, Phillips & 

Wilson (1980) over the fitting function for a^j^ provided 

by Verner et al. (1996): 



n^ycaf/ = r(Mg°) = 6.8 x lO-^^^ s'^ . (A4) 

Although the radiation field for the gas associated with a 
GRB host galaxy will likely have a different spectral index 
than that of the Milky Way, the ionization rate is domi- 
nated by the photon density at /li^ w 8 cV and, in any case, 
uncertainties in the atomic data limit the accuracy of this 
analysis to factors of a few. Rearranging Equations A1-A4 
to express Ue in terms of the observed iV(Mg°)/A''(Mg+) 
ratio we reproduce Equation 4. 

One can derive a similar expression for Fe. In this case, 
radiation recombination dominates for T < 8000K, 

o^'^^iT) = 8.945 X 10-9 [(Vr/0.04184(VT/0.04184+ 1)°-^«^^(1 + ^/T/KSblTW^y-^^^^ 

(A5) 

For the ionization rate we perform the same calculation as 
above (Equation A4) and find 

r(Fe+) = 1.7 X 10-i°(G/Go) . (A6) 
Finally, here are the relevant expressions for Carbon: 

a^+(T) = 7.651 X lO'^ [ (a/T/0.001193( ^2^/0.001193 + l)0 i983(i _^_ ^T/9.334 x I0i2)i-8027 

(A7) 

and 



r(C+) = 6.6 X 10-^°(G/Go) . (A8) 
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